Spin blockade in a charge-switchable molecular magnet 
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We consider the effect of adding electrons to a single molecule on its magnetic properties and the 
resulting transport fingerprints. We analyze a generic model for a metal-organic complex consisting 
of orbitals with different Coulomb repulsions. We find that by modulating the charge of the molecule 
by a single electron the total spin can be switched from zero to the maximal value supported by the 
added electrons, S = 3/2. The Nagaoka mechanism is responsible for this charge-sensitivity of the 
molecular spin. It is shown that fingerprints of these maximal spin states, either as groundstates 
or low-lying excitations, can be experimentally observed in current-spectroscopy. as either spin 
blockade at low bias voltage or negative differential conductance and complete current suppression 
at finite bias. 

PACS numbers: 



•4-* 



C 

O 

o 



o 

O 

in 
o 
-i— > 



T3 

c 

o 
o 



X 



Introduction — Recent experiments on metal- 
organic grid complexes, consisting of rationally designed 
ligands and metal ions as building units have exhibited 
interesting electrochemical Q, Q and magnetic 0, 0, Q 
properties. By self-assembly the metal ions and ligands 
arrange in a rigid, highly symmetric grid. Due to their 
different nature, electron orbitals can often be roughly 
attributed either to the metal-ions or the ligands. 
Typically the orbitals on the (organic) ligands have 
7r symmetry whereas the d metal orbitals split in an 
e.g. octahedral ligand field into subshells with local 7r 
and a symmetry. In the case of a fully occupied d v 
subshell tunneling between ligand -k orbitals and metal 
d(j orbitals is weak due to their different symmetry. 
Such a separation into metal ion and ligand units has 
been used successfully to describe the low-temperature 
intramolecular spin coupling of Co- [2 x 2] grids 0, Q 
and Mn-[3 x 3] grids [j, |fj for a fixed charge state as 
well as the electrochemical properties of (Mn, Fe, Co, 
Zn)-[2 x 2] 8] and Mn-[3 x 3] grids 0. For poly-pyridine 
complexes it is well-known 9] that ligands as well as 
metal ions can be reduced. Which type is preferred 
depends on chemical details which can be controlled, 
mainly by substitution of metal ions and changing 
the ligand. Here we analyze a phenomenological low 
temperature model for a [2 x 2] grid molecule consisting 
of four ions and four ligands, Fig. ^ For this particular 
structure we show that (i) the molecular spin can be 
switched by the charge and (ii) the spin-splitting appears 
in tunnelin g sp ectroscopy. The well known Nagaoka 
mechanism [lOj becomes effective for certain numbers 
of added electrons. For strong onsite interaction the 
delocalization of an extra hole/electron relative to 
half-filling (favoring a fully polarized background of all 
other electrons) dominates over an antiferromagnetic 
superexchange. In the context of band-magnetism the 
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FIG. 1: Grid molecule: small/large circles represent metal 
ion/ligand orbitals. 



relevance of this mechanism is limited due to its lattice 
type dependence and its strong charge sensitivity. 
Only for a single additional electron or hole relative 
to a half-filled band, the spin-polarization effect can 
be guaranteed. In small single-molecule devices these 
obstacles can be overcome. Firstly, the advanced ratio- 
nal design of supramolecular structures allows complex 
"lattice" types to be realized 0, 0. Secondly, due to 
the Coulomb blockade effect one can modulate the total 
charge of a molecule by a single electron 0, EH E3- 
In the case of sufficiently strong short-range interaction 
on the ligands (relative to the ligand-ion tunneling) the 
ground state spin is maximally increased from S — to 
the maximal value supported by 3 (or 5) extra electrons 
on the molecule, S = 3/2. The magnetic properties of 
the molecule can thus be controlled electrically 0, . 
In single electron tunneling transport this leads to spin 
blockade at low bias voltage. Even for a low-lying max- 
imal spin excitation negative differential conductance 
(NDC) effects and complete current suppression at finite 
bias voltage occur. Similar models with two types of 
electron orbitals have been studied for the description 
of the neutral-ionic transition [Tsl ] in organic crystals, 
in the context of ferroelectrics and superconductivity in 
transition-metal oxides Jljj and recently in the context 
of exotic Kondo effects due to dyna mical symmetries in 
multi quantum-dot systems |201 121| . 

Model — We consider a model of a grid-complex with 
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four metal and four ligand sites and one orbital per site 
(Fig. (IJ. For simplicity, we assume that the metal ion 
can only be occupied virtually. The strong ligand field 
separates electron accepting d orbitals energetically from 
ligand orbitals. Additionally, Coulomb repulsion on the 
ions is typically much stronger than on the ligands. The 
following Hamiltonian captures the features of the elec- 
tronic degrees of freedom: 

H mo \ = Hj + H\_ + Hm, (1) 

h t = j2J2 tA l a ^ + h - c - ( 2 ) 

4 

2 

+ W ^ n 3 n J+2 (3) 

i=i 

4 

H M = ^{ENi + UN^Nn) (4) 
i=i 

Operators and variables (except t) in lower/upper case 
relate to the ligands/metal ion and all indices run 
from 1 to 4 cyclically. denotes a summation over 

neighboring metal ions i and ligands j. The operator 
a] CT creates an electron on ligand site j with spin er, 
iij a = a^cijo- and nj = J2a n j<?- Similar definitions hold 
for the metal ions: A ia ,N ia = J^ la A la ,~N l = J2 a N i<?- 
The tunneling term ((2J) describes hopping between 
ligand and metal ions and is assumed to be independent 
of i and j due to molecular symmetry. The ligand-part 
of the Hamiltonian © consists of an orbital with 
energy e, the Coulomb repulsion terms on the ligand 
(it) and between adjacent (i>) and opposite ligands 
(w). Due to decreasing overlap with distance we have 
u > v > w. Hamiltonian Q describes the isolated 
metal ion orbitals with energy E. Here we only consider 
the short-range interaction U because the d orbital 
overlap between two ions is typically much smaller 
than that between two ligand orbitals. In Fig. these 
interactions are schematically indicated. We study the 
parameter regime where the first eight extra electrons 
occupy four equivalent ligand centered orbitals. Such 
a sequence has been well-documented for a number of 
grid-molecules 8]. In our model this is the case when 
the charge excitations of the ligand lie below the ones of 
the metal ion: e < e + u < E < E + U. The metal ligand 
charge-transfer barrier A = e — E suppresses the direct 
hopping of extra electrons from the ligands to unoccu- 
pied metal ions: |A| 3> t. The fluctuations of the metal 
orbital occupation around zero can be treated using a 
Schrieffer- Wolff transformation [2^]. We obtain an ex- 
tended Hubbard model iJ e fF = J2(jk) J2a *eff a\ a a ka + H L 
on four ligand sites with an effective hopping matrix 
element i e fj = — 

Addition energies and spin states — We first analyze the 
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FIG. 2: Splitting AEn ag for n — 3 as function of u for v = 
2.25t, w = It, A = -lOt, tai = 0.051 For u t h w 4.15* = 83i e ff 
the ground state has maximal spin. Inset: ground state spin 
as function of the number of electrons n added to the ligands 
for u < u t h (solid black line) and it > u t h (dashed blue line). 



electronic spectrum of H e ff for different charge sectors. 
In the case u > v > w/2 > £ e ff the first electron reduces 
one of the four ligands, say j ' = 1. The next one occupies 
the opposite ligand j = 3 in order to minimize the 
Coulomb interaction. The next electrons reduce the 
adjacent ligands j — 2,4. This sequence is repeated 
for the next four electrons, each time doubly occupying 
a ligand orbital. The gaps in the addition spectrum 
w, 2v — w,w,u,w,2v — w,w (extra energy required for 
the next electron) thus directly relate to geometrical 
features of the grid ||. These electrostatic parameters 
correspond directly to the size of the Coulomb diamonds 
obtained in transport experiments. Now we discuss the 
ground state spin as successive electrons are added to 
the molecule. Due to superexchange processes electron 
spins on neighboring ligands tend to couple antiferro- 
magnetically. This leads to an alternating sequence of 
S = and S = 1/2 as electrons are added (inset Fig. EJ. 
However, for sufficiently large u > w t h and fixed t, resp. 
t e ft (Fig. |2| the ground state spin for n = 3, 5 is enhanced 
from the noninteracting value S = 1/2 to the maximal 
possible value S — 3/2. Because double occupation 
is suppressed, a single hole/electron (relative to the 
half-filled state n = 4) can maximally gain kinetic energy 
when the background of the other electrons is fully spin 
polarized. This is the underlying mechanism for the 
Nagaoka theorem [To| . The ferromagnetic alignment 
(due to complete delocalization) competes with the 
antiferromagnetic spin coupling (due to superexchange 
processes between neighboring occupied sites). Which 
process dominates depends on the strength of the onsite 
repulsion u relative to fixed hopping. The interactions 
v, w tend to increase the threshold value w t h for fixed 
teff p3j . The gap between the Nagaoka state and the 
lowest excited state saturates at AE^ag ~ 2t e ff (Fig. 0) 
independent of u due to the kinetic origin of the effect. 
In order to attain an observable effect one should thus 
have sufficiently large values of both u and t e ff. For 
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FIG. 3: dI/dV(V g , V) grayscale plot (white/black ^ 0) for 
u = 5t,v = 2.25i, w = t, A = -lOt, T = 2 10 _4 i. 



example, w can be increased by a chemical modification 
of the ligands which draws charge density into the 
ligand LUMO orbitals. Taking typical parameters Q 
|A| w leV,i w 10 _1 eV, AE Nag « lCr 2 eV, we estimate 
u t h ~ leV which is reasonable. 

Transport — The charge sensitivity of the total spin 
is observable in the single electron tunneling-current 
through the molecule. To demonstrate this, we consider 
the Hamiltonian H = H res + H e ff + H mo \- res , employ- 
ing units h = e = ks = I- The electrodes r — L,R 
are described as electron reservoirs with electrochemi- 
cal potentials fi r = fj, ± V/2 and a constant density 
of states p: H res = J^kar e ka r cl ar c kr7r . The tunnel- 
ing term ff mo |_ re5 = {^) 1/2 T, k ajr *i c Lr a ^ + h - c - de- 
scribes charge transfer between electrode and molecule 
(symmetric tunneling barriers) . T is the overall coupling 
strength between leads and the molecule and defines the 
current scale. We assume that tunneling is only possi- 
ble through two "contact" ligands, namely t\ = = 1, 
otherwise 0. We have checked that this choice does not 
cause effects due to orbital symmetry as discussed in [24| 
by trying also j = 1,2. The coupling to a gate elec- 
trode is included in a shift of the single particle energies 
e — ► e — aV g . In the weak tunneling regime (T <C T) 
the effect of the leads can be incorporated in the transi- 
tion rates S S)S ' = J2r ^st' + ^t>7> between the molecular 
many-body states s,s': 

= Tj2tf(Es-E s ,)\Y^tr(s\al\s')\ 2 

" 3 

= r^/-(E s - J B s oi^*i(«K<rk'>l 2 - (5) 

Here /+ is the Fermi function of reservoir r and /~ 
1 — /+. Importantly, the matrix elements include the 
calculated many-body wavefunction of the molecule and 
the spin selection rules. From the stationary master 
equation ^2 s i{^s,s'Ps' ~ S S ', S P S ) = we obtain the non- 
equilibrium occupations P s of the molecular states s and 
the resulting stationary current which may be calculated 
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FIG. 4: Same parameters as Fig. [3] except u — 4.05t and 
around a different charge degeneracy point. 

at either electrode r = L,R: 

i r = -J2(K:t'Ps'-Kr s p s ). (6) 

Due to the presence of a maximal spin state, either as 
ground or excited state, spin blockade and NDC effects 
occur, respectively po| . 

Maximal spin ground state — For u > u t b and fixed 
t the ground state for n = 3, 5 is a Nagaoka state 
(Fig. El- Transport involving groundstates n <-» n + 1 
for n = 2, • • ■ , 5 is completely blocked for small bias 
and low temperature (T,V < AE^ 3g ): the Coulomb 
diamonds in a differential conductance versus (V g ,V) 
plot do not close (Fig.|2J). Since the ground state spin is 
either (n = 2,4, 6) or 3/2 (n = 3, 5) the tunneling rates 
between neighboring ground states vanish: even when 
ground state transitions are energetically allowed, due to 
the spin selection rule AS = 1/2 transport is completely 
blocked in the weak tunneling limit. However, when 
temperature or voltage are increased, such that the first 
excited state with appropriate spin can be accessed, 
current begins to flow. 

Maximal-spin excited state — Depending on the gate 
voltage, NDC and even complete current blocking can 
occur at finite bias voltage |25| when the Nagaoka state 
is the lowest spin-excitation for n = 3, 5 (i.e. u < u t h, 
FigEJ. The typical result is depicted in Fig. 0] Near 
the charge 3 h 4 degeneracy point (and 4 <-> 5) two 
NDC lines (black) with negative slope appear which 
are due to the low lying 5 = 3/2 Nagaoka state. We 
first discuss the lower NDC line using the left scheme in 
Fig. At the charge degeneracy point (V g /t ~ 6.58 in 
Fig. the current sets on because a transport channel 
is opened, namely n = 3,S = 1/2 (2-fold degener- 
ate) <-> n — 4, S — 0. Increasing the bias voltage by A4 
results in a gain in population of the Nagaoka state via 
the n = 4, S = 1 excited state. Since the Nagaoka state 
cannot decay (strictly) to the n = 4, S = ground state 
the number of transport channels is therefore decreased 
from two to one leading to the lower NDC effect. 
Further away from the degeneracy point {Vg/t < 6.56 
in Fig. 2J, the lower NDC line turns into a conductance 
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FIG. 5: Minimal set of states for the lower NDC effect (left) 
and the complete current suppression (right) in Fig. 0] 



peak and simultaneously the ground-state transition 
line below it disappears. This is due to a complete 
population inversion between the ground and excited 
state for n = 3 which already occurs inside the Coulomb 
diamond. When the transition from n = 4, S = 
to the third excited state n = 3,S — 1/2 lies in the 
bias window, the Nagaoka state is occupied starting 
from the ground-state n = 3, S — 1/2 via the cascade 
of single-electron tunneling processes indicated in the 
right panel of Fig. [3] It is fully occupied because the 
escape rate from n — 3, S — 3/2 relative to that from 
n = 3, S = 1/2 ground-state is suppressed by a factor 
~ e -(A4-A 3 )/T < 1 (Fig. 03 . The upper NDC line with 
negative slope in Fig. 0] is caused by the occupation of 
the high lying maximal-spin state n — 4, S = 2 which 
cannot decay to states with one electron less and higher 
spin [23. This state can already be reached at low 



voltages only due the presence of the low-lying Nagaoka 
state S = 3/2 at n = 3. 



Conclusion — We have shown that in particular types 
of single-molecule devices the Nagaoka spin-polarization 
mechanism is relevant since the Coulomb blockade effect 
allows the controlled addition of single electrons (in 
contrast to the case of bulk magnets). The magnetic 
properties of the molecule may thus be switched by a 
gate voltage. The fingerprints of this charge-selective 
stabilization of maximal spin states (cither a ground or 
a low-lying excited state) are observable in the tunneling 
current. We also have investigated an extension of the 
model considered here with a spin degree of freedom 
added to each of the four metal ions |2£|. In this more 
complicated case the cooperative spin-polarization effect 
(of both Nagaoka-origin and direct ion-ligand exchange) 
compete with the antiferromagnetic superexchangc 
(induced by ion-ligand hopping). It is found that the 
Nagaoka state is also relevant here and that its transport 
fingerprints are the same as demonstrated in this work 
for a generic model. 
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